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Abstract 

The microlensing objects, Machos, recently observed in the halo of our Galaxy, can be 
interpreted as dense neutralino objects, neutralino stars, produced by the gravitational in- 
stability of neutralino gas. Taking the mass and radius of these objects from microlensing 
observations we calculated the diffuse gamma-ray flux produced in neutralino-neutralino 
annihilation inside the objects. The resulting flux is many orders of magnitude higher 
than the observed one. 
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The recent observations of microlensing objects in the halo of our Galaxy [[TJ show that 
the Macho mass is O.46to;5yM at 68% CL. This large mass of lenses raises the question of 
how they can escape from the observations by the Hubble telescope. The fraction of the 
total halo mass in the form of Machos, according to the observations ffl], is £ = 0.50lo'|o- 
If this fraction is higher than 50% a new question arises as to why the fraction of non- 
baryonic cold dark matter, which is needed for the large scale structure formation, is so 
small in our Galaxy 

In some recent papers a very interesting idea about the nature of these microlensing 
objects was put forward. These objects can be dark matter formations around a singu- 
larity, produced during the non-linear stage of the fluctuation growth (|| and references 
therein). The authors call these objects neutralino stars (NS). If this idea is correct, the 
microlensing phenomenon is a natural result of fluctuation growth in the neutralino gas, 
and the aforementioned possible difficulties connected with the baryonic nature of Machos 
disappear. 

In this note we shall study the gamma-ray production in the neutralino stars in an 
approach more general than the model [[|. Namely, we shall adopt for the neutralino 
stars the mass, the radius and their space density in the halo being restricted by the 
microlensing observations. For the density distribution inside NS we shall use p oc r~ L8 , 
which is a basic theoretical result for this kind of objects |§. We shall consider the mass 
of neutralino as a free parameter, though in the work |2| the neutralino mass is restricted 
as 5 GeV < m x < 50 GeV. 

The diffuse gamma-ray flux from neutralino stars in the halo is 

J 7 = ^—risRhN^ , (1) 

47T 

where n s is the number density of neutralino stars in the halo, Rh is the radius of the 
halo and iV 7 is the production rate of photons in the NS. 

The density of NS in the halo can be estimated by assuming that the fraction £ of the 
total mass of the halo, Mh, is in the form of NS: 

where M x is the mass of the NS. 

Let us now estimate the production rate of photons iV 7 in an NS. 

The density of neutralinos inside the NS depends on the distance as r~ 18 || and can 
be expressed through the mass M x and radius R x of the NS as 



with a negligible correction due to a modification of the r 18 dependence at very small 
distances. The gamma-ray production rate is 



iV 7 = 2r]^Aix(ov) j n 2 x (r)r 2 dr (4) 

R, 
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where (av) is the cross-section of neutralino-neutralino annihilation and r) n o(m x ) is the 
multiplicity of neutral pion production in neutralino-neutralino annihilation. 

We have assumed in (4) that the neutralino star has a core of radius R c . Inside the 
core the neutralino density is constant or rising more slowly than r _L8 , and we neglected 
the gamma-ray production there, which is 5 times smaller than outside it for the case of 
constant density. 

Using Eqs. (1), (3) and (4) we obtain the diffuse flux as 

0.9 £M h M x /R x \ 06 (av) 
^"47^° ml \Rj R x Rl' [b) 

To proceed further we should evaluate two parameters: core radius R c and annihilation 
cross-section (av). 

Let us start with the annihilation cross-section (av). This quantity is of relevance both 
for the determination of the flux (see Eq. (5)) and for the evaluation of the neutralino 
relic abundance. In the usual expression (av) = a + bx, where x = m x /T, valid in the 
non-relativistic limit, a takes contribution from the s-wave only, whereas b contains both 
s and p-wave contributions. The relic abundance is given by 
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n x h 2 = 8.67 ■ lO- 11 ^ , Ge ^ /0 ^ 2 ( G ) 

x n) 1 ax f + (l/2)bx) 

where Nf ~ 100 is the number of degrees of freedom at the epoch of the neutralino decou- 
pling and xj = m x /Tf ~ 1/20, Tf being the temperature at the decoupling. Situations 
in which this approximation is not valid are discussed, for instance, in ||. 

Let us note that the annihilation cross section in Eq. (5) has to be evaluated at 
the present temperature (To) and then for xq — 10~ 7 . This implies that 7 7 critically 
depends on the size of the a term in (av); any suppression effect in the a term entails a 
large depletion in the diffuse flux emitted by NS. As will be shown later, the flux 7 7 has 
the tendency to largely (by many orders of magnitude) exceed the experimental bound. 
Thus in what follows we adopt the conservative approach of analysing with special care 
situations in which the term a is significantly suppressed. 

The evaluation of (av) requires the inclusion of the full set of annihilation final states 
(// pairs, gauge-boson pairs, Higgs-boson pairs and Higgs gauge-boson pairs), as well as 
the complete set of Born diagrams (Z, Higgs, squark, neutralino and chargino exchanges). 
Furthermore one should include loop contributions, leading to gluon final states, which 
may be relevant in annihilation processes at the present temperature. Analytical expres- 
sions and numerical codes for the annihilation cross section are available (for a complete 
set of references, see [Q). 

As discussed above, the case favourable for neutralino stars as microlensing objects 
is the one in which the s-wave annihilation, described by a, is strongly suppressed in 
comparison with the p-wave contribution given by b. This occurs in the gaugino and 
higgsino dominated states when the neutralino is lighter than the IV-boson and the final 
states are fermionic pairs. In the case of strong s-wave suppression (bxf ^> a) one obtains 
from (6): 
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b~ 1.4- l(T 35 -^7 cm 2 . (7) 

Our way of writing the numerical factors in the previous expression is motivated by 
the observation |7| that for all cosmologically successful models the cold dark matter 
contribution to Q is given by QcDMh 2 = 0.2 ± 0.1. 

Turning now to the calculation of (erf) at the present time, we shall evaluate a using 
b as given by Eq. (7) and the ratio a/b, which is almost free from particle-physics uncer- 
tainties following the procedure suggested in Ref. |§. Let us first consider the case when 
mj < m x < mw, where m& is the mass of the 6-quark. The dominant contribution to the 
cross section is given in this case by bb final states. 

For the gaugino-like case the ratio of a/b is given |J by 

a/6=6 . 3 . 1(r *(W^r) 2 //( ,), (8) 

where mioo = m x /100 GeV. The function f(y), where y characterizes the neutralino 
composition, is given by || 

, _ 567 - 10% + 1242?/ 2 - 12y 3 + 2023?/ 4 

and is bounded as 

8 < f{y) < 120. (10) 
For a higgsino-like neutralino, the a/b ratio is ||: 

4.i • io- 6 ( mb/5 GeV V. (ii) 

o V m ioo / 

Notice that the smallness of the a/b ratio in both the gaugino- and higgsino-like neutralinos 
is due to the well known suppression effect in the s-wave annihilation channel (proportional 
to the square of the mass of the final-state fermion) due to helicity properties [|J. 

We recall that for the case of a neutralino of a mixed gaugino-higgsino composition 
the Higgs exchange in the annihilation cross section is usually important and the size of 
the a/b ratio depends critically on the values of the Higgs masses. For instance, for a light 
A Higgs boson (Ma < 2m x ) the a term dominates, at decoupling, over the b term in the 
annihilation cross section. This is due to the fact that, for the Higgs exchange diagrams, 
also the 6-term is proportional to m 2 because of the Yukawa couplings. The case of s-wave 
suppression here can occur only at the pole Ma = 2m x (see Ref. 0), but in this case 
Q x h 2 is less than cosmologically needed. 

Therefore, to minimize the gamma-ray production in the NS we consider the Higgsino- 
like neutralino with the a/b ratio given by Eq. (|TT|) . Then using Eq. (7) we obtain, for 
the cross-section: 

(av) w 1.7 • 10- 30 m7 2 cmV 1 , (12) 
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Note that this is the smallest annihilation cross-section we can obtain. Since b is fixed by 
£l x h 2 (Eq. (7)), the value of the cross-section in Eq. (|i~2| ) is determined by the suppression 
of the s-wave contribution. 

For a neutralino heavier than the W^-boson, the situation is different. Among many 
channels with weak gauge bosons and Higgses in the final states, there are those where 
the s-wave is not suppressed at all and those where it is strictly forbidden. Thus the 
average suppression is not strong. 

A discussion of the core radius R c is in order now. We shall consider first the case of a 
pure neutralino star, with no baryon contamination. The core radius can be estimated by 
equating the rate of neutralino accumulation in the core, 4irR 2 n c u r , and that of neutralino 
annihilation there, (4ir/3)Rln 2 (crv ), where n c is the density of neutralinos in the core and 
u r is the bulk streaming velocity of neutralinos towards the centre (we shall omit the 
index r in the following discussion): 

R, = . (13) 

n c {av) 

To estimate u, let us consider the Euler and Poisson equations ||, which determine the 
flow of the neutralino gas: 

du du d(p 
dt dr dr 

A<p = AnGp, (15) 

where <fi is the gravitational potential, p is the gas density and G is the gravitational 
constant. 

The solution (3) corresponds to the quasi-stationary regime, when du/dt can be ne- 
glected. Putting (3) into Eq. ([15]) one finds 



90 _ GM X f_r\ °"' 
dr~ R\ \R X ) 



(16) 



Now it is easy to solve Eq. (|14"D , which gives 



u(r) = v / 5(^) ai « // , (17) 
R x 



where Uff is the free-fall velocity for an NS: 



I2GM X 
R x 



Another way to estimate u(r) is just to assume that it is everywhere the free-fall 
velocity relative to the mass inside the radius r. In this case we obtain 

(19) 
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to be compared with Eq. (|Hj). 

Now by inserting Eq. (|17D into Eq. (|13D one finds 

R c = 28R^ m5 m^ 33 M°f 5 cm, (20) 

where Ru is the radius of the neutralino star R x in units of 10 14 cm and Moi is M x in 
units of O.IMq. 

We shall further minimize the flux (||) by using the following set of parameters : 
m x = mw = 80 GeV, 77^0 = 10, taken from e + e~ — > hadrons data at ^/s = 80 GeV, 
R max = 200 kpc and consequently M h = 1O 12 M , £ = 0.1, M™ in = M macho = O.1M , 

R max = 1 . 1Q 14 ^ = 2 Q . IQ-30 cm 3 g -l &nd R c = \ . 1Q 2 C m. 



The diffuse flux at £ 7 > 70 MeV is 



jm*n _ x . 10 7( 10 2 cm /R c f- b cm-'s-^r" 1 , (21) 

almost 12 orders of magnitude higher than allowed by observations of SAS 2 and EGRET 
(J 7 < 2 • 10~ 5 cm~ 2 s -1 sr -1 at £7 7 > 70 MeV [|TIJ). As was mentioned above the value of 7 7 
in Eq. (|21|) is already a conservative estimate in terms of neutralino composition. Even 
adopting the very extreme phenomenological assumption that in the annihilation at the 
present time the b term is dominant (i.e. a/b < 10 -7 ), the value of the gamma flux would 
only be reduced by a factor ~ 40, as compared with the value in Eq. (21), and then it 
would still be largely above the experimental bound. 

Why is the flux in Eq. (|21~D so large? It is because we compressed the considerable 
fraction of the neutralino mass in the halo in the dense NS, taking the radius of neutralino 
star R x < 10 14 cm from the focusing condition ||. It is instructive to compare the flux 
(5) with the case when neutralinos are distributed homogeneously in the halo. The diffuse 
gamma-ray flux is 

A° m = ^/V*ft> ( 22 ) 
where the rate of gamma-ray production /3 7 is 

/3 7 = (n;° m ) 2 2 V (^> 

and n^ om is the homogeneous density of neutralinos in the halo. The ratio of fluxes from 
NS and from a homogeneous distribution of neutralinos is 

This ratio can be easily rewritten as 

rNS „ „ ^ / d x 0.6 



^7 < Px > NS ( R\ 



(24) 



lit (Px) hom \ Rc 

where one can easily recognize the compression factor discussed above as the ratio of av- 
erage neutralino density (p x ) in the NS and in the halo. Numerically this ratio is ~ 10 17 . 
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The flux given by Eq. ([H]) is obtained for a case of a pure NS. In reality the neutralino 
star should unavoidably have some baryonic contamination, because dissipative baryonic 
gas is streaming to gravitational potential minimum created by neutralino gas and is 
accumulated there. This phenomenon was clearly recognized and considered by Gurevich 
et al. (see and references therein). We would like to note here that heating of the 
baryonic gas due to neutralino-neutralino annihilation (which was not taken into account) 
is very important for the dynamics of the baryonic accretion in vicinity of a singularity. 
For the mass of the baryonic object M b the ratio M b /M x > 1CT 3 was suggested in Ref. 0. 
We shall use M b = 0.01M X pT| . As a simple estimate shows, for this mass the neutralino 
trajectories are modified by the baryonic core at distances R c < 0.03R X . However, it 



is questionable how the space density of neutralinos is changed. In Refs. [|T^, it is 
shown that the density distribution of a non-dissipative gas with baryonic contamination 
is the same as for a non-dissipative gas. Even assuming that the core radius is very large, 



R c = 0.1R X UTTH, the calculated flux turns out to be P ~ 3 cm s sr , still 5 orders of 
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magnitude larger than allowed by observations. 

In conclusion, our analysis shows that the interpretation of the observed Machos or 
even of a small fraction of them as neutralino stars results in too large a diffuse gamma- 
ray flux, many orders of magnitudes higher than the observed one []. 
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